Nonrelativistic Fermions in Magnetic Fields: a Quantum Field Theory Approach 
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The statistical mechanics of nonrelativistic fermions in a 
constant magnetic field is considered from the quantum field 
theory point of view. The fermionic determinant is computed 
using a general procedure that contains all possible regular- 
izations. The nonrelativistic grand-potential can be expressed 
in terms polylogarithm functions, whereas the partition func- 
tion in 2-1-1 dimensions and vanishing chemical potential can 
be compactly written in terms of the Dedekind eta function. 
The strong and weak magnetic fields limits are easily stud- 
ied in the latter case by using the duality properties of the 
Dedekind function. 



I. INTRODUCTION 

The motion of fermions in magnetic fields is an old 
problem of quantum mechanics that plays a role in a 
wide variety of physical problems, ranging from neutron 
stars Q to the quantum Hall effect Q . From the point of 
view of quantum field theory, several studies have been 
performed, related to the chiral anomaly Q, effective ac- 
tions 1^-^ and applications to anyons systems 0. 

In this paper we would like to study another aspect of 
this problem, namely, following the approach proposed in 
Q , we shall discuss the statistical mechanics of nonrela- 
tivistic fermions embedded in a constant magnetic field. 

Our approach requires the explicit computation of a 
fermionic determinant, and this is carried out in terms 
of polylogarithm and Dedekind functions. Additionally 
our calculation, in a particular case, allows us to study 
the strong and weak magnetic field limits quite easily, 
due the modular invariance of the Dedekind function. 
The paper is organized as follows: in section 2 we review 
the approach proposed in in section 3 we formulate 
the problem in the context of the path integral method, 
in section 4 we present an explicit computation of the 
nonrelativistic grand potential, and in section 5 we study 
the partition function in the case of 2-1-1 dimensions. Our 
conclusions are given in section 6. 
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II. EFFECTIVE ACTION AT LOW ENERGIES 

In reference a method was proposed to compute 
effective actions at low energies, based in a path inte- 
gral derivation of the Foldy-Wouthuysen transformation. 
This method can be used for nonrelativistic fermions as 
well. Let us consider the lagrangian 



where Ip = P 

Following the procedures developed in 
the origin of the energy by the rescaling 

to write the lagrangian as 

£ = 0(i^-m(l-7o))0. 



(1) 

we redefine 
(2) 

(3) 



Now we decompose the spinor into 'large' {lp) and 
'small' (x) components, in whose terms 

L = ip^ iDaip + [iDQ + 2m] x + « (i • Dx 

+X^iS- D^, (4) 

where the Dirac representation for the 7-matrices has 
been used. 

The next step is to diagonalize the lagrangian in (Q) 
by means of the change of variables 



X =x+[«^o + 2m] ^ id -Dip, 

ipUa ■T>[iDo + 2m]~'^. 



(5) 



This change of variables has a Jacobian equal to unity 
and the effective lagrangian, under this transformations, 
becomes (omitting the primes) 

C = (p'^[iDa + a ■'D{iDo + 2m)~'^ a -D] p, 

+ X^[«^o + 2m]x. (6) 

This lagrangian describes the non-local dynamics of 
the fermions in terms of two-components spinors. One 
should note that p and x decouple and, after expanding 
(iI?o + 2m) ^ in powers of l/m, the partition function 
becomes 



Z = j Vx^Vxe'^^ j Vip'fVipe'^^, 



(7) 



1 



where 



S^^ I d*x^^ [ iDo + -^D2 + -^a • B] ^ + 0(1/7 
J '■ 2m 2m ■' 



(8) 



Sip is the nonrelativistic action for fermions interacting 
with a magnetic field. S^^ is the action associated to the 
lower contribution of the spinor and can be neglected in 
the nonrelativistic limit. 



III. NONRELATIVISTIC FERMIONS IN 
MAGNETIC FIELDS 

Using (0) and (||) one can explicitly study the nonrel- 
ativistic quantum field theory of fermions at finite tem- 
perature. 

The partition function associated to this problem is 



^ dct [iDa + T^T)'^ + ■ B + fi], (9) 
2m 2m 



(10) 



where g stands for the electric charge of the fermions, fj, 
is the chemical potential, and the infinite product runs 
over all the eigenvalues A„ of the operator that appears 
in @, 

[iDo + ^T)^ + ^<J-B + fi]q^n^XnK, (11) 
2m 2m 

subject to the usual antiperiodic boundary conditions in 
the imaginary time direction. 

For a constant magnetic field one can choose the gauge 

A = (-Boy, 0,0), 
Ao = 0, 

in which case the fermionic determinant can be computed 
from 



(12) 



and thus the partition function (^ becomes 

= det[idt - H+ + fi]dct[idt - R- + ^i], (13) 



n n 



where the can be read off from (p2|). 

Each determinant in ( p^ ) is evaluated by explicitly 
solving the eigenvalue equation (^2|) by means of the 
Ansatz (/'^(x,t) = /±(x)T(i), which yields 



it{t) - Q.T{t) = 0, 



(15) 
(16) 



since the operator 2^D^ ± -^Bq is time-independent. 
Equation (15) is just Schrodinger's equation for the Lan- 
dau problem, whose eigenvalues are known, 

Et = n-\t = {n + ^±^)u; + ^-^^, (17) 

with uj — gBo/m. The equation for T(t) has a solution 
only if 



n^n,^^ —{2m+l), 



(18) 



where T is the period and m an integer, in virtue of the 
antiperiodic boundary condition on T{t). 
Thus, the eigenvalues in (p^) are given by 



A 



r 



(2m -M)-^;,^ 



(19) 



The statistical mechanics of the fermionic system un- 
der consideration is described by the grand potential - 
basically the logarithm of the partition function after go- 
ing to Euclidean space, i.e. replacing T by i(3, where 
(3 — 1/T is the inverse temperature. The logarithm of 
the partition function is 



oc oc 



^ogZ^ = J2 / dPz[logiX+^J + log(A„^„)], 

71—0 m— — oo 
oc „ 

/ dp,[i+ + L,;], (20) 

n=0 

where are infinite sums defined as 

oo 



) 



(21) 
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Although the series in (^l|) are divergent, they can be 
computed by using a definite regularization prescription. 
Let us start considering the divergent series 



am + b). 



(22) 



m— — oo 



where a and b are constants. 

The second derivative of L{a, b) is 



(PL{a,b) 
d62 



E 



1 



(am -I- bY ' 
TT^ csc^(67r/a) 



m— — DO 
^2 „„„2 



(23) 



so that upon integration one finds 
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L(a, b) = log ( e'=i+^"^ si 



(24) 



with 



where ci and C2 are two arbitrary integration constants 
In our case, a and b can be read from ( |l9|) and (21) 
Therefore 



= log 



-i^t] cosh (^/3:|^ 



zy + (-c± + |)£;±+log(l 



(25) 

(26) 



Clearly, the constants cf and parametrize the arbi- 
trariness of the regularization. Since the Euclidean space 
effective action must be real, we choose ™ such a way 
that (26) has no imaginary part. 

At this point we can compare our result (^6|) with the 
nonrelativistic limit of the general result given in Q for 
the effective action. The contribution 

cf - + {-4 + l)Et 

indeed corresponds - after adding the analogous positron 
contribution- to the first term, Tr|5|, in equation (18) of 
reference and the dependence on the arbitrary con- 
stants reflects the regularization that needs to be 
performed in order to define the trace. For instance the 
choice ct = = is consistent with ^-function regular- 
ization Q. 

The second contribution, log(H-e^'''^'> ), which is finite 
and independent of the constants 0^*^21 coincides (up to 
a factor f3) with the grand potential in nonrelativistic 
statistical mechanics. 



The function S{A, b) has the following Taylor series in 
the variable A: 



5(A,6) = log(l + A) + ^ J J_ 4", 

n— 1 ^ ^ 

nil - e-"'') 

n— 1 ^ ^ 



(32) 
(33) 



Since has the form a"e "'•'^^ , the pz integral can be 
performed explicitly to yield 



(1 - e-"'') 



(34) 



Expanding now 

(1 _ e-"'')-! = ^ (e-"'')^ 

we find 



k=0 



r dpzS{A^,b) = -.I^Y.Uy2{~a^e-''), (35) 

^ ^ k=Q 

where we have introduced the polylogarithm function, 
defined as the analytic continuation to the whole complex 
z-plane of the series 



IV. GRAND-POTENTIAL 



(36) 



In this section we compute the nonrelativistic grand 
potential {Q) directly from its thermodynamical defini- 
tion, rather than deriving it form the partition function 
(Z). 

The nonrelativistic grand potential is given by 
^{P,f^) = -^ E / dpzJ2\ogil + e-''<), (27) 

' 1=+.--'-°° n=0 

= n+i(3,fi) + n-i/3,f,), (28) 

where g is the degeneracy factor g — [gB^)/ Ati'^)V . 
Clearly, we need to consider the following object 



defined for Ul < 1 and s G C. Therefore, 



S{A,h) = Y,\og{l + Ae 



n=0 

In terms of S we have 



i7±(r,^) = / dpzS{A^,puj) 



(29) 



(30) 



n+{T, fi) = V—-^ ^ Li3/2(- 



4^2 ^3/2 



.g/3Mg-(fe+l)/3c^)^ 



/c=0 



(37) 



and 



n-{T,^,) = y|-^^^Li3/2(-e''^) + J7+(T,^). (38) 

It is instructive to check that our result for Q, = 
ri+ -I- reduces to the correct result in the limit of 
zero temperature (/? — > oo). Depending on the sign of 
fi — kbj, the argument of the polylog function will tend 
either to zero (in which case Li3/2(0) = and there is 
no contribution) or infinity. To get the corresponding 
contribution in the latter case we need the asymptotic 
behavior of Li3/2(z) as z ^ oo. This can be obtained 
from Joncquiere's relation 119], 



3 



Li,(z) + e'^" Li,(l/z) 



log 2; 



271-i 



In particular, for a; > we have 



Li3/2(-e'^") 



3i7r/4 



r(3/2) 

2 /33/2 



11 (3x , 
2' 2 2^' 



3r(3/2)' 



p3/2 



(39) 

(40) 
(41) 



3^ 



because of the asymptotic behavior C(— i,g) - 
of the Hurwitz zeta function. With r(3/2) = \/tt/2, we 
obtain the correct nonrelativistic result 



37r2 ^ 



2m 



^ (M - fcc)^/^ - 



3/2 



fc=0 



(42) 



Here [xj is the floor of the real number x. 

A similar calculation can be performed in 2+1- 
dimensions, but then the result one obtains is nothing 
more than the starting expression (p8|), with the corre- 
sponding 2-|-l-degeneracy factor ^2 = {9Bo/2tt)L'^ , 
being the area of the 2-dimensional quantization box: 



172+1 (r, m) 



92 



fc=0 



log(l 



1 + e/3Mg-(fc+i)/3'^ 



(43) 



V. PARTITION FUNCTION IN 
2+1-DIMENSIONS 

In section we compute the partition function in 24-1- 
dimensions. In order to do that, one can proceed as fol- 
lows: the product dl3) is first written as 



Zip — {ZoY 



(44) 



in view of the degeneracy 32 of each Landau eigen- 
value. Using ^-function regularization the reduced parti- 
tion function Zq can be put into the form 



cosh 



n=0 





cosh 




2 




2 



(45) 



Using the explicit form of the eigenvalues for the Lan- 
dau problem in 2-f 1 dimensions (pz = 0) from (17), ( p5| ) 
can be written as 



^0 = n cosh(6'n + S+) cosh(6lr 



(46) 



where 



Defining q 

becomes 



e = 








s- = 






2 ' 


, z = 





(47) 



n 



n=0 



-,-20 



2710 S-^+S- 



(1 + Z<Z")(1 + Z(Z") 



(48) 



This is most general expression for the (reduced) par- 
tition function in 2+1-dimensions. 

A particular situation arises when the chemical poten- 
tial is zero. In such a case (Esf) becomes 



Z. 



r]{2id/Tr) 
rj{i6 /-k) 



where 



(49) 



(50) 



is the Dedekind eta function ||T^ , with q = e" 

It is not difficult to see that formula (|4^) reproduces 
result ( p3[ ) in the case of vanishing chemical potential. 
Indeed, from the definition of the eta function and writing 
1 - = (1 _ (,")(! + g») we find 



^0 = 9^/^^11(1 + '?")'' 

n=l 

which, with q = e~^^ leads to 

^ 00 
logZ^ = --(3uj + 2 ^ log (1 + e-"''") 



(51) 



(52) 



This coincides with the result for —l3il2+i from (|42_ 
at /i = 0, up to terms due to the regularization used. 
The high temperature limit (or equivalently the weak 
magnetic field limit) corresponds to q 1. where the 
Dedekind function can be approximated by |13| 



1/2, 



(53) 



and, therefore, the partition function in the high temper- 
ature limit becomes 



^0 



cosh^ 



(54) 



n=0 



Moreover, as the Dedekind function satisfies the mod- 
ular property 



4 



1/2 

V 



(55) 



one can compute the partition function in the strong 
magnetic field hmit (or low temperature limit) directly 
from (|5^). Indeed, by using (j5^) and ( p5[ ) we find 



7^ 



cosh ( — — 
\l3uj 



(56) 



VI. CONCLUSIONS 

In this paper we have studied the motion of nonrela- 
tivistic fermions in a constant magnetic field, following a 
quantum field theory approach. The main point is the 
calculation of the fermionic determinant by using a gen- 
eral method that contains all the possible regularizations 
of the infinite product, in the form of particular choices 
of the coefficients cf2- 

In addition, we obtained explicit expressions for the 
nonrelativistic grand-potential, in terms of a series of 
polylogarithms, and for the partition function in 2-f 1 di- 
mensions and zero chemical potential, in terms of the 
Dedekind eta function. In the latter case, the modular 
properties of the Dedekind function allow us to relate 
the strong and weak magnetic field limits. This relation 
could be useful in the dynamics of the quantum Hall ef- 
fect. 
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